Production of scalar particles in expanding spacetime 
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In this paper, we investigate cosmological particle production using quantum field theory (QFT). 
We will consider how production of scalar particles can occur in an expanding universe. By in- 
troducing a time-dependent energy parameter representing the time evolution of the universe, the 
initial vacuum state will be excited. Consequently, creation of particles is present. Here, our focus 
is mainly creation of minimally coupled scalar particles in Minkowski spacetime. 
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I. INTRODUCTION 



Particle creation is an important topic in modern 
physics. The concept was given birth after the discovery 
of the inconsistency of vacuum. One of the first investi- 
gations of particle production by an expanding spacetime 
was first carried out by Bernard and Duncan in 1977. It 
showed that the concept of vacuum was incorrect, since a 
vacuum state does produce particles in a spacetime that 
undergoes some form of expansion. The creation of min- 
imally coupled scalar particles, or free bosons, will be 
considered in this paper. There are many important rea- 
sons to consider scalar particles. These particles are cer- 
tainly generated only by the geometry of the spacetime. 
In the calculation of the number of these peculiar bosons, 
we will see that these bosons mostly are zero-mode par- 
ticles. Somehow, this may indicate the so-called spin-2 
gravitons. The idea may also lead to an explanation of 
the vacuum energy divergences of bosons and fermions. 
Clearly, particle creation occurs in pairs, i.e. particles 
and anti-particles will be created as a pair. From the 
calculations in this paper, we will see that the initial 
vacuum state can be represented by a 2n-excited final 
vacuum state. 

In Birrell 01 , this topic has been investigated in a two- 
dimensional Friedman-Robertson- Walker spacetime with 
a suitable line metric. In his analyse, the expansion of 
the spacetime is represented by the conformal scale fac- 
tor. In this paper, however, we will investigate the same 
problem in term of QFT. In the derivation, produced par- 
ticles will be identified as excitations of an initial vacuum 
state of a simple harmonic oscillator. Considering Klein- 
Gordon fields in an expanding spacetime represented by 
a time dependent energy term, minimally coupled scalar 
and neutral pairs of particles will be produced. Certainly, 
the evolution of the field is given by the so-called Bogoli- 
ubov transformation and the number of produced parti- 
cles is represented by the absolute square of the Bogoli- 
ubov coefficient, |/3kP- 



II. THE MODEL 

Here, we consider a simple model to see how excita- 
tions of a vacuum state can occur in a usual Minkowskian 
spacetime at distant past and distant future. Actually, 
since the term vacuum is not adequate, it should be re- 
placed by ground state. But for simplicity, we will not 
be concerned by this, so in this model we may consider 
the ground state as a vacuum state. The ground state 
at t ^ — oo is associated with an initial vacuum state 
which evolves in time. Indeed, the simplest way to start 
the investigation of this topic is a simple massless har- 
monic oscillator. For massless particles, there is equality 
between the momentum and the energy. Later, in consid- 
eration of massive bosons, the energy uj{t) can be written 
as a function of momentum k and mass m in a quan- 
tum theoretical way. In our case, with the assumption 
of a time dependent energy the corresponding Hamilto- 
nian acquires the same time dependence. By finding the 
probability of the transistion from the initial to the final 
vacuum state, excitations of the initial vacuum state will 
clearly be present. Obviously, this will also appear in 
calculation of the equation of the motion. Excitations of 
the initial vacuum state will appear in the calculations, 
which simply can represent production of particles due 
to time evolution. 

Our starting point is a massless harmonic oscillator 
with a time dependent Hamiltonian given by 



Hit) = \p[tf 



(1) 



where q{t) is the scalsr field and p{t) is the momentum 
of the field. Both depend on the time t and no other 
parameters in this simple model. The frequency is chosen 
to be 



Aw^ tanh(7i). 



(2) 
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which is our ansatz for the time evoulution of this model. 
This represents an expanding spacetime where far past 
and future are expressed as asymptotic regions t — > — cx) 
and t — > -|-c». Since we are working with a massless 
field and the ansatz of the frequency is not a function 
of the spatial coordinate x, spatial translation invariance 
is still valid and is a symmetry of this spacetime. It is 



2 



thus convenient to consider the field as a quantum field, 
yielding 



q{t) =au{t) + a^u*{t). 



(3) 



Here, a and are the annihilation and creation operators 
and u{t) is the time dependent function of the field. The 
Langrangian is L — /2 — LiP{t)q^ /2 and the canonical 
momentum is p = dL/dq = q. Here, we have set the mass 
parameter in the common harmonic oscillator equal to 1 , 
and 7 is a parameter adjusting the shape of the energy 
evolution. For the very early universe, i.e. t — > — oo, we 
let uj{t) uji, and for the distant future, uj{t) —* 102, 
i.e. t — > 00. To simplify analysis we only consider the 
asymptotic cases of i — s- ±00, corresponding to the case 
of sudden limit, i.e. 7 ^ 00. Hence, the frequency will 
change its form at the point of i = 0, i.e. 



_ \ujI - Au;^ = LUi for t <0, 
^ \Lj'^ + Auj^ = UJ2 for t > 0. 



(4) 



From the Euler-Lagrange equations, we obtain the fol- 
lowing equations 



q{t)+^{t)\{t)=0, 



(5) 



where uj{t) is given by Eq. The derivative with 

respect to t is denoted by an overdot. We consider q{t) 
given by Eq. ||3J) to be a time dependent quantum field 
on the form 

q{t) = ^={a,e-'^'* + a\e'^''), t<0 

where and a are the creation and annihilation oper- 
ators and the prefactor comes from normalization. The 
field must obey the conditions of continuity q{t = = 
q{t = 0+) and q{t = 0") = q{t = 0+). The conditions 
give the following equations. 



and 



UJl 



ai + a\ 



:(ai - a\) 



a2 + 



W2 



;(a2 - 4) 



\J2lo\ \j2Zj\ 
respectively. We define a new parameter r\ such that 



(6) 



(7) 



(8) 



Writing Eqs. © and iQ in matrix form gives 



1 1 



1 



1 



02 



(9) 



By defining 



1 1 



the creation and annihilation operators for t > from 
Eq. @ can be writing as 



02 




Eq. can be writting in a more compact form, reading 



0-2 



cosh r\ sinh r\ 
sinh r] cosh r\ 



ai 



(11) 



where 



coshr/ = 
sinh 77 = 
The inverse form of Eq. H1U|) is 




cosh 77 — sinh t] 
— sinh ?; cosh ?; 



0-2 



(12) 



Clearly, we have obtained an operator relation between 
the early universe and the present universe. The intro- 
duction of cosh 77 and sinhry above can be justified by 
associating with the relations 



t A ^ 

X a' and — a. 

ox 



(13) 



For an arbitary vacuum state the following relations 
are prevailing 



ai|r!>i -0, 
a2\n)2 = 0. 



(14) 



Obviously, Eq. (|12|l gives ^cosh 7702 — sinh 7703^ = 

0. Using the associations p3f) . the annihilation operator 
02 can be represented hy d/dx and the creation operator 
a\ by X in quantum mechanics. Fortunately, the initial 
vacuum state can be expressed as an excited state of 
reading 



^— tanh 7] 



\n)2 



(15) 
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III. THE VACUUM STATES 

The conditions of the definition of a vacuum state for 
t ~* ±00 are given in Eq. (^H). From Eqs. and (fT^ . 
it is possible to find the relation between the vacuum 
states. From Eq. H12() . we have the relation 



(cosh 77 a2 — sinh?7 a2)|ri)i = ai|r2) 



0. 



(16) 



Obviously, a2\^)i is non-zero. According to Eq. lfTC|) we 
can make an ansatz of the initial vacuum state, reading 



(17) 



where M is the normalization constant. Using the associ- 
ation Eq. H13(l . the normalization constant can easily be 
find. Utilizing series expansion of the exponential func- 
tion in Eq. (|17|l . it can be written as a squeezed vaccum 
state with ^ = — tanh?/ in the form, 



\o = {i-\e\)"'Yl 



(-O'^TM! 



fc=0 



|2fc), 



where \2k) = 
state H17II reads 

00 



|ri). Accordingly, the initial vacuum 



— tanh 77 
2 ' 



2n 



'2n\ 



n=0 



■ tanh rj 



2n 



I2n\ 



where we have used the coefficient of probability 



/2nJ 



C2n = 



— tanh 77 
2 ' 



(18) 



Clearly, the initial vacuum state |r2)i is a kind of a 27i 
excitated state of |ri)2, since a'2(^-^^2 ~ |277)2. Utilizing 
the Stirling's formula 



(19) 



can simplify the coefficient C2n in a more satisfying form 
V2^ = (47r77)i/42"n"e-". Hence, = \/2^77"e-". 
The prefactor in front of the tanh-function in Eq. I|18l) 

/2^ 1 _ (47r)i/" 



can be written in the form, — j— 7^ — , — , 
We can use the condition of normalization. 



for n — !■ 00. 



-1/2 



(—4)". Consequently, 



E 



2n \ f tanh 



E 



77 / y 2 

-1/2" 

77 

1 



(— 2tanh77)'^ 



z , (21) 
^ 1 — tanh rj ^ 

The normalization constant can be found by using Eq. 
(Pn|l and the relation 7V"^ = l/i/l - tanh^ 7/. Hence, 
the normalization constant of the initial vacuum state 
reads M = I/v^cosEt/. The initial vacuum state can thus 
be expressed by the present vacuum state, yielding 



1 



\/cosh rj 



From the definition of every excited state, \n) = 
-^(a^)"|r2), the vacuum state |ri)i is clearly a excited 

state of 1^2)2, 



^cosh ?7 



E 



/2n} f\ 

— f-tanhr,) |2r7)2, (22) 



where |277)2 = ((4)^"/%/2^) |^^>2• 

According to Eq. Hll() . we have the relation 

(cosh7/ai + sinhr/a{)|f7)2 — a2\^2 — 0. (23) 
which implicitly gives the commutation rule reading 



02,6 



^ tanh r\ {cii^)^ 



_ „(itanh7,)"(at)2 



(24) 



Obviously, Eq. is a consequence of Eq. H23|l . since 
this yields 



^ tanh J] (02)^ 



] = E^(^tanh7y)"[«2,(4n 

°° 1 1 

= - (tanh 7y) (- tanh 7y) '" 

X 7^4(4)^" 

(777 -1- 1)! 

1 

777!'2 



tanh 7; 4 E ;;^(|tanh77)"(4; 

m— 



2m 



by using the commutation rules [A, BC] ~ \A,B\C 
B{A, C] = and [02, = M {AY^^^- 



(f]|r!)i = l. 



(20) 



to rewrite Eq. p7|) into a geometric series. The 

/ 277. 

rewritting is clear by utilizing the identity 



IV. THE PROBABILITY OF EXCITATIONS OF 
THE VACUUM STATE 

It is interesting to find the probability of excitations 
in the initial vacuum state at t — > — 00 due to the time 



4 



evolution. In our case the excitations represent the num- 
ber of produced particles in the harmonic oscillator, and 
the definition of the initial vacuum state, is the 

state at t — *■ — oo. Vacuum state here means the low- 
est energy state of the Hamiltonian. The transition from 
to 1 17) 2 will initiate particle production depending 
on the energy relation, r/ defined in Eq. (|HJ|. According 
to Eq. H22|l . there is evident that 2n quanta will be cre- 
ated during the transition. Hence, the probability of the 
transition is found by placing the dual vector 2{^\ into 
Eq. (|22|l . Due to the normalization (|2Uf) . the probability 
of the transition is found to be 



P2„{\n)i ^ \2n)2) 
= \2{2n\n),\^- 



coshr; 



2n 
n 



(tanhry)^"^ 



(25) 




where n = 0,1,2,.. and 



2n 



(n!)2 



is a kind of a 



binomial coefficient, i.e. the number of ways of picking 
n unordered outcomes from 2n possibilities. We see that 
P2n has a maximum value of rj for each n. If Aw = 0, 
i.e. r] = 0, Po is maximal at 77 = 0. However, if we 
consider two particles are produced by the vacuum state 
when t 00, P2 is maximal at 7? ~ 1.146, i.e. the 
present vacuum energy must be 1.318 greater than the 
vacuum energy at the early universe. 

By considering 77 and Eq. H25|l . the vacuum state in 
the early universe, |fi)i, can be identified as an excited 
state of the present vacuum state, \^)2- 



« 



1=0 



Tl = 1.146 



1.444 



It = 1.628 



FIG. 1: The 77-value for 2n particles are produced by 
when t — > 00. If there is no energy difference ij] = 0) between 
these two vacua, and r2)i, the Hamiltonian cannot be 

time dependent and therefore no particles can be produced. 
However, if 77 > 0, there will produced particles in pairs. 



V. THE GENERAL SOLUTION OF THE 
EQUATION OF MOTION 

In order to find the expression of the field qit)., we must 
solve the Euler-Langrange equation |SJ 



q{t) + {ujI + tanh7i) q{t) = 



(26) 



FIG. 2: The probability for zero, two, four, eight particles 
are produced by the time evolution of the state t —* 00. 

When 1 7; I increases, the probability of more produced particles 
becomes greater. 



with 



lim q{t) = — == 



Then, the question is how the field behaves when the 
time evolves to infinity, i.e. q(t —^ 00). Let us generalize 
the field. We start with writing the field for t —00 in 
usual quantum field theoretical form, reading 



q{t) = aiui{t) + a|u*(t). 



(27) 



where ui{t) is a time dependent orthonormal basis. We 
want to find ui(t) which satisfies the differential equation 
Eq. (|^. By inserting Eq. I^Tj) into we certainly 

obtain 



ili{t) + (cjo + Atj^tanh7i)ui(t) = 0. 



(28) 



In order to identify Eq. H28|) with an equation with known 
solution, we introduce the new variable 



Tit) = i(l+tanh(7i)). 



(29) 



Thus, the double derivative in time can be written in the 
form 



= 47^ 



r{l-T)]^-^+T{l-T){l-2T)] — 



dr 



Writing Eq. H28|l with this new parameter, one certainly 
recognizes that this is a hypergeometric equation which 
reads 



\ [r{l - t)]2— + r(l - t)(1 - 2r)]— \ ui(r 



-t- WgUi - Aa;^(l - 2t)ui{t) = 0. 



dr 



(30) 
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Clearly, r goes to when i — > — cx) and r goes to 1 when 
t ^ oo. The solution of Eq. (|30|l has a typical form 

oo 

reading mi(t) = t"^ ^ it„T", where is a parameter 

independent on r. The hypergeometric equation H30|l 
has three singular points t — 0, t ~ 1 and t ~ oo. We 
define the two time independent aymptotic vaulcs of Lu{t), 



reading ujf = ujq 



Aw^ and lo^ — cj. 



for t ^ — c 



and i — > oo respectively. The index equation for r = 
by assuming ui ~ t'^° at the singularity t = reads 



(31) 



The solution acquires the form ui ^ r^*(~) for r — > 
0. To obtain the standard form of QFT, we choose the 
negative sign of the indices of the solution. Hence, 



Ml 



(32) 



For r = 1, the solution of the field will go like ui 
(1 — ■ The index equation thus becomes 



vi{vi - 1) + ^1 + 



w 



0. 



(33) 



This gives us the solution of Eq. (|30|) for the case of 
distant future, reading 



Ui 



(1-r) 



(34) 



Analog to l|32() . we choose the right sign of the index of 
(|34|l that satisfies the standard form of quantum field the- 
ory. Since (1 — ^ exp(— 27t), it is therfore naturally 
to choose the positive sign of the index of (|34|l . Hence, 
the solution can be written in the form 



ii = (l-T)'(w)r-<^)/(T), 



(35) 



where /(r) is again a hypergoemetric function, /(r) can 
be solved by inserting it back to the orginal hypergeo- 
metric differential equation H3U|) . Inserting (|35|l into Eq. 
()30|l certainly gives 



1 



r(l 



1 - 



IUJ2 

1 



-{LUl 

7 



f' + Tf^O, 



(36) 



where ' denotes the derivative with respect to r. T is 
a r-independent term coming from the last term in Eq. 
(|30|l and by deriving of /(r), yielding 



T 



27 



47 



Certainly, Eq. H3t)|) is a hypergeometric differential equa- 
tion with the form 

T{l~T)f" + [c-{a + b+ 1)t] f'~abf^O. (37) 



By comparing the general form of hypergeometric equa- 
tion (|37|) with our equation (|36|l . we simply obtain 



27 



i{uji + UJ2) 



27 



1 + I'D ~ i'l , 



7 



= 1 



(38) 

(39) 
(40) 



These expressions enter in the hypergeometric function 
/(r) — const F{a,b : c,t) 0. Clearly, the solution of 
the hypergeometric differential equation contains singu- 
larities in T = and t — 1, yielding ui{t) ~ t^^=^^)(1 — 

r)(l^)F{a,b;c;T). 

For T ^ 0, i.e. t —^ —00, we obtain the same form we 
had to begin with Eq. since F{a,b;c;T) goes to 1. 
The frequency mode becomes 



ui{t) ~ e 
where we have used 



-iuJi t 



t^-oo (41) 
(42) 

Considering the singular point t — 1, i.e. t — *■ 00, the 
frequency mode ui{t) will obtain a time evolution given 
by a hypergeometric function. 



where Fum is given by 



Fiini = lim F{a, b; c; r) 

a:— *1 + 



lim 1 



(43) 



— lim 



r(c)r(c - a - b) 



r(c- a)r(c- 6) 
r(c)r(a + 6-c) ,^ 



i^(a, b;a + b - 



r(a)r(6) 

xF{c — a.c — b\c 



- a- 6+ 1: 1 



r)] 



1; 1-r) 



(44) 



This function becomes unity when r — > 0. The relations 
which enter in Eq. I|44|l are 



c — a — 5 = 21^2/7, c — a = 1 
i(w2-t^i) 



J(W2 - UJl) 



b-c=- 



27 

icJ2 



27 7 
Hence, for the case of t — > 00, the solution will go like 



(45) 



(1-r) 



-27* 



(46) 



when we transform t back to t. Here, the hypergeometric 
functions F{a, 6;a + 5 — c— 1; 1 — t) and F{c — a,c — 
b; c ^ a — b + 1; 1 — t) will go to 1. Thus we obtain two 
terms in the solution of ui for t —>■ 00, containing time 
independent coefficients, K_ and K^, in each of them. 
Certainly, we get the forms of solution for those two cases 
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t ^ zLoo. By using Eqs. (|41(l and 1)46(1 . the frequency 
mode acquire the forms, Ui{t) = Ce^''^^* for t — > —oo 
and ui{t) = Ci e''^^^K+ + C2 e-"^^^K^ for t +00. We 
have used the relations given in Eq. (|45|) and defined the 
two time independent products of gamma functions as 



and 



K- = 



V{c)V{c- a-b) 
T{c-a)T{c-h) 



T{c)T{a + h- c) 



(47) 



(48) 



The constants C, Ci and C2 in the frequency must have 
a form so that the field q satisfies the communtation re- 
lations, [a, flt] — 1 and [a, a] — [a^,a^] = 0. However, the 
constants also must be in such way that the Lorcntz in- 
variance of the field is conserved in the continuous form; 
e.g. if the frequencies loi and uj2 are dependent on the 
wave vector k, the integral of the field over k must be 
Lorentz invariant. Consequently, the constants thus ob- 
tain the forms C = l/\/2wi, Ci ^ C2 ~ l/\/2tj2, which 
satisfy the hypergeometric equation Ij^UI) . 

Now, we will write down the asymptotic solutions of 
the field. For t — > ~oo, we have that r ~ e^''*. Thus the 
frequency mode has the form given by H41I) is obtained, 
since 



T 27 



-iuJi t 



And for t 
(1 - r)^^ 



00, (1 - r) 



-27* 



It thus gives 



which is consistent with H43|l . Con- 
sequently, the solution of the field will obtain the usual 



form for t - 
form for t 



ui{t) 



where J^K 



—00, i.e. q{t) 
00. Indeed, 



oiui +a[ul and a modified 



1 g-ii^lt 



/2wl 
1 



-00. 



(49) 



and 



+ are the so-called Bogoli- 

ubov coefficients for the massless bosons, denoting by a 
and /3 respectively. From H49|l we see that the expression 
(t) obey the Bogoliubov transformation with the form 



Ui = au2 + (3u 



2; 



a U2 + p U2, 



(50) 



where m — 



_1 ^—icJit 

— ~ 1 



i ~ 1,2. The annihilation opera- 



tor will do the same with the form 



* /3* t 

a 02 — p flj. 



Pa2. (51) 



In order to consider massive bosons, we include mass 
term in our ansatz of the evolution of the energy uj(t), 
yielding 



(52) 



where i = 1,2. As for the massless bosons, this ansatz 
of the frequency is not a function of the spatial coor- 
dinate X. Thus, spatial translation invariance is still 
a symmetry here. The solution (|49|l will be the same, 
but now Lu{t) and the solution of the field ui{t) for the 
massless case include a k-dependence, i.e. uj{]i,t) and 
w(k, i) — (l/v^2tJi) exp (— iwii) exp (ikx) . 

Thus, the Bogoliubov coefficients for the massive min- 
imal coupled bosons read 



r(i 



7 ' 



iji r(- 



27 



)r(i 



27 - 



K- (53) 



^1 ^(^^^)^(l + ^ 27 



(54) 



By using T{x) r(l — x) = sin-wx 3' ^^^^ obtain the 
moduli of the Bogoliubov coefficients 



a- 



sinh^(7rti;+ /j) 



sinh(7r(jJi/7) sinh(7ra;2/7) ' 



sinh^ (ttoj- /7) 
ih(7ra;i / 7) sinh(7rijj2 / 7) ' 



(55) 



(56) 



where we have defined the energy differences, oj- = 
i(w2 — oji) and uj+ = \{u)2 + wi)- Since the commu- 
tation relation of the operators for t — > cx) is [ai,a|] — 
(|ap — |/3p)[a2,a2] = 1, the normalization condition of 
the Bologoliubov coefficients reads 



1. 



(57) 



Grafically, the coefficients can be represented by Fig|31as 
a function of momentum |k|. Certainly, the Bogoliubov 
coefficients we have found for a harmonic oscillator with 
the ansatz of the energy (|52|l are exactly the same results 
in Ref. |3|]. Birrell (1982) has given a derivation in parti- 
cle production by considering this topic as a cosmological 
problem describing by a metric for the Robert- Walker- 
universe with a conformal scale factor which corresponds 
the time evolution of our model of harmonic oscillator. 



VI. THE PROPERTIES OF THE BOGOLIUBOV 
COEFFICIENTS 

The Bogoliubov coefficients at and /3k that we have 
found in Eqs. 1)53(1 and 1(54(1 can be evaluated as 



Otij = [Uii, U2j) 



(58) 



where x is the four-vector. The Bogoliubov transforma- 
tions l(5U() and l(51() are both diagonal and isotropic. Their 
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modes of the (vacuum) field in that way we can identify 
excitations of the initial vacuum state as created 

particles. 



VII. ADIABATIC EXPANSION 



o.a- 



0,6- 



0,4 



0,2 



Q l I 



|k| 



FIG. 3: The graphical representation of the Bogoliubov co- 
efficients for the produced bosons. In this figure we compare 
the numerical coefficients with the expected analytical val- 
ues, for a particular choice of parameters mi = f, 7 = 10 and 
7712 = 10. 



coefHcients prossess the following properties 

a-ika^k - Ak/3j*k = %7 (59) 

k 

a^k/^jk - AkCKjk = 0. (60) 



It follows that two spaces based on the two choices 
of modes ui and U2 are different as long as /3k 7^ 0. 
As we have considered ealier, a vacuum state will 
not be annihilated by the annihilation operator 02, i.e. 
a2|ri)i = /3k|l)i 7^ 0. The state |l)i is a excited state 
of the initial vacuum state Indeed, by using the 

Bogoliubov transformations the expectation value of the 
number operator N2 — a\a2 of U2-frequency mode parti- 
cles in the initial state reads 



v{n\N2\n), 



(61) 



which states that the initial vacuum state with frequency 
modes U2 has been excited with |/3kP excitations. Simply 
stated, we can say that |/3p given in Eq. (|55|l is the 
number of particles produced in the ii2-mode due to the 
time evolution oiLo{t). It can thus interpret as production 
of minimally coupled scalar particles, i.e. bosons, in a 
time dependent harmonic oscillator. According to Eq. 
(|61|l . the absolute square of /3k represents the expectation 
value of the number of particles present in the final state. 
From another point of view, the time evolution that we 
have assumed in our model can be considered in a pure 
cosmological context. Hence, the evolution represents the 
smoothly expanding motion of an aymptotically static 
universe. Indeed, the expansion excites the frequency 



Considering an adiabatic expanding spacetime, the dif- 
ference of the energies uj- goes to zero. The number 
of produced particles represented by Eq. H56|l will thus 
vanish. Indeed, no particle can be created. This is a 
conformally trivial situation, i.e. a conformally invari- 
ant field propagating in a spacetime that is conformal 
to Minkowski spacetime. Hence, the particle creation oc- 
curs only when the conformal symmetry is broken by the 
presence of mass providing a length scale for the theory. 
The particle production can be identified the coupling of 
the spacetime expansion to the quantum field through 
the mass. The expansion regarded as a changing gravi- 
tational field will be an energy supplier to the perturbed 
scalar field modes. When we consider creation of field 
quanta with the changing gravitational field, it is natural 
to describe the particles as being produced in the asym- 
potically static regions. Certainly, particle production 
can not occur in the asymptotically static regions. How- 
ever, if represents a final density of particles which 
were not present in the state we would except that a 
measurment taking place at an intermediate time during 
the expansion era, would unveil a particle density in the 
interval from zero to the number of particles is produced 
given by the absolute square of /3k. 

This do not live up to scrutiny, however, if the space- 
time is curved, there will no good definition of particles is 
generally available. Because of the special symmestry of 
the FRW spacetime, we can identify comoving observers 
who see the universe isotropically expands. In this case, 
it is reasonable to consider produced particles in the ex- 
pansion region as the excitation of comoving particle ob- 
server. 

Let us get back on track and arrive at the mathemati- 
cal description. The adiabatic case certainly gives limits 
of Eqs. (|17|l and regarded to the Bogoliubov coeffi- 
cients (|53|) and (|54|l . i.e. 



K_ 



r(i 



^)r(^) 



^^2) 



r(i) 

r(i-^)r(-i^) 
r(-f )r(i-i^) 



27 



(62) 
(63) 



since 0^2 = t^^i = in ths case. In an adiabatic ecpanding 
spacetime, it is appropriate to assume no energy change, 
. i.e. LUi — L02. The field q(f) is approximately equal 
for both asymptotical cases. According to Eqs. H62|l and 
(inSJ, the solution (gnj acquires the form 



ux{t) 



(64) 
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This clearly states no change in of the field. The field evo- 
lution satifes the normalization condition (|57|l according 
to Eqs. ^ and yielding P - 0^ = 1. The number 
of created particles, . i.e. number of excitations, is rep- 
resented by 1/3 p. In an adiabatic expansion, the creation 
rate approachs to zero when the expansion rate go to 
zero. Here, it is appropriate to consider 7 as the expan- 
sion rate since it is the parameter of the time evolution. 
When 7 — > 0, Eq. H56|l goes thus like 



— -27rwi/7 



0, (65) 



which states that the number of created quasi-particles 
rapidly approaches zero in an adiabatically expanding 
spacetime. 



where we have assumed that the energies for both asymp- 
totic cases are positive. Since the field expansion changes 
in the case, the number of created particles can not be 
neglected. For sudden limit, the expansion is able to 
generate particle production. However, the bosons which 
are been created due to a sudden expansion of the space- 
time is nevertheless very small. According to |/3kP, there 
will be produced mostly zero-mode particles, i.e. k = 0. 
Therefore, in any calculation of the bosonic vacuum en- 
ergy, which is a sum over k, the created bosons will thus 
be invisibles in the final {i.e. t ~f 00) vacuum energy 
calculation, since the number of created bosons is mainly 
zero-mode particles. 



VIII. SUDDEN LIMIT 



IX. CONCLUDING REMARKS 



A. Sudden limit 

For sudden limit the expansion rate increases dramat- 
ically in the time evolution, i.e. 7 — > 00. Certainly, the 
time parameter T(t) given by (|29|l approaches to a step 
function. The terms on the hypergeometric function, 
and , will go like 



7 -i{uji - UJ2) 
iuj2 
1 



27 
'i{uji 



:(1 



OJ2 ' 



-IUJ2 



27 



1 . 1-^1 s 



(66) 
(67) 



In the sudden limit approximation, things become more 
dramatically. The initial vacuum state will accordingly 
be excited (or de-excited) due to the enormous change 
of the spacetime. The initial state energy uji evolves to 
bj2 for t 00. According to the Eqs. H66|l and (|67|l . the 
solution (|49|l acquires the form 



t 



Ul{t) 



Ml 



-)e 



2V- ■ ui2' 
iuj2t 



-00 



00 



(68) 



Fortunately, the field evolution satifes the normalization 
condition (|S7I) . reading 



In this papir, we have given a derivation of particle 
production in an expanding spacetime by using a simple 
harmonic oscillator model and the Bogoliubov formalism. 
We have shown that particle creation in a time-depending 
harmonic oscillator can occur in a spacetime with a suit- 
able ansatz of the time evolution, corresponding to the 
motion of the expansion of the universe. Since the very 
beginning of the universe underwent an enormous expan- 
sion, this spontaneous production of particles played a 
vigorous part at that time. According to the calculation, 
created particles are equivalent to excitations of the ini- 
tial vacuum state. In this model, we have only considered 
bosonic case of the cosmological particle production. 
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